Biological control networks suggest the use of biomimetic sets for combinatorial therapies

Cells are regulated by networks of controllers having many targets, and targets affected by many controllers, but these "many-to-many" combinatorial control systems are poorly understood. Here we analyze distinct cellular networks (transcription factors, microRNAs, and protein kinases) and a drug-target network. Certain network properties seem universal across systems and species, suggesting the existence of common control strategies in biology. The number of controllers is ~8% of targets and the density of links is 2.5% \pm 1.2%. Links per node are predominantly exponentially distributed, implying conservation of the average, which we explain using a mathematical model of robustness in control networks. These findings suggest that optimal pharmacological strategies may benefit from a similar, many-to-many combinatorial structure, and molecular tools are available to test this approach.Comment: 33 page

Tubular secretion of creatinine and kidney function: an observational study.

BackgroundPrior papers have been inconsistent regarding how much creatinine clearance (CrCl) overestimates glomerular filtration rate (GFR). A recent cross-sectional study suggested that measurement error alone could entirely account for the longstanding observation that CrCl/GFR ratio is larger when GFR is lower among patients with chronic kidney disease (CKD); but there have been no validation of this in other cohorts.MethodsTo fill these gaps in knowledge regarding the relation between CrCl and GFR, we conducted cross-sectional and longitudinal analysis of the Modification of Diet in Renal Disease study (MDRD) and African American Study of Kidney Disease and Hypertension (AASK); and cross-sectional analysis of a clinical dataset from the Mayo Clinic of four different patient populations (CKD patients, kidney transplant recipients, post kidney donation subgroup and potential kidney donors). In the cross-sectional analyses (MDRD, AASK and Mayo Clinic cohort), we examined the relation between the CrCl/iothalamate GFR (iGFR) ratio at different categories of iGFR or different levels of CrCl. In the MDRD and AASK longitudinal analyses, we studied how the CrCl/iGFR ratio changed with those who had improvement in iGFR (CrCl) over time versus those who had worsening of iGFR (CrCl) over time.ResultsObserved CrCl/iGFR ratios were generally on the lower end of the range reported in the literature for CKD (median 1.24 in MDRD, 1.13 in AASK and 1.25 in Mayo Clinic cohort). Among CKD patients in whom CrCl and iGFR were measured using different timed urine collections, CrCl/iGFR ratio were higher with lower iGFR categories but lower with lower CrCl categories. However, among CKD patients in whom CrCl and iGFR were measured using the same timed urine collections (which reduces dis-concordant measurement error), CrCl/iGFR ratio were higher with both lower iGFR categories and lower CrCl categories.ConclusionsThese data refute the recent suggestion that measurement error alone could entirely account for the longstanding observation that CrCl/GFR ratio increases as GFR decreases in CKD patients. They also highlight the lack of certainty in our knowledge with regard to how much CrCl actually overestimates GFR

A Stochastic Multiscale Model of Cardiac Thin Filament Activation Using Brownian-Langevin Dynamics.

We use Brownian-Langevin dynamics principles to derive a coarse-graining multiscale myofilament model that can describe the thin-filament activation process during contraction. The model links atomistic molecular simulations of protein-protein interactions in the thin-filament regulatory unit to sarcomere-level activation dynamics. We first calculate the molecular interaction energy between tropomyosin and actin surface using Brownian dynamics simulations. This energy profile is then generalized to account for the observed tropomyosin transitions between its regulatory stable states. The generalized energy landscape then served as a basis for developing a filament-scale model using Langevin dynamics. This integrated analysis, spanning molecular to thin-filament scales, is capable of tracking the events of the tropomyosin conformational changes as it moves over the actin surface. The tropomyosin coil with flexible overlap regions between adjacent tropomyosins is represented in the model as a system of coupled stochastic ordinary differential equations. The proposed multiscale approach provides a more detailed molecular connection between tropomyosin dynamics, the trompomyosin-actin interaction-energy landscape, and the generated force by the sarcomere

Mechanisms of transmurally varying myocyte electromechanics in an integrated computational model

The mechanical properties of myocardium vary across the transmural aspect of the left ventricular wall. Some of these functional heterogeneities may be related to differences in excitation–contraction coupling characteristics that have been observed in cells isolated from the epicardial, mid-myocardial and endocardial regions of the left ventricle of many species, including canine. Integrative models of coupled myocyte electromechanics are reviewed and used here to investigate sources of heterogeneous electromechanical behaviour in these cells. The simulations (i) illustrate a previously unrecognized role of the transient outward potassium current in mechanical function and (ii) suggest that there may also exist additional heterogeneities affecting crossbridge cycling rates in cells from different transmural regions

Modeling effects of L-type ca(2+) current and na(+)-ca(2+) exchanger on ca(2+) trigger flux in rabbit myocytes with realistic T-tubule geometries.

The transverse tubular system of rabbit ventricular myocytes consists of cell membrane invaginations (t-tubules) that are essential for efficient cardiac excitation-contraction coupling. In this study, we investigate how t-tubule micro-anatomy, L-type Ca(2+) channel (LCC) clustering, and allosteric activation of Na(+)/Ca(2+) exchanger by L-type Ca(2+) current affects intracellular Ca(2+) dynamics. Our model includes a realistic 3D geometry of a single t-tubule and its surrounding half-sarcomeres for rabbit ventricular myocytes. The effects of spatially distributed membrane ion-transporters (LCC, Na(+)/Ca(2+) exchanger, sarcolemmal Ca(2+) pump, and sarcolemmal Ca(2+) leak), and stationary and mobile Ca(2+) buffers (troponin C, ATP, calmodulin, and Fluo-3) are also considered. We used a coupled reaction-diffusion system to describe the spatio-temporal concentration profiles of free and buffered intracellular Ca(2+). We obtained parameters from voltage-clamp protocols of L-type Ca(2+) current and line-scan recordings of Ca(2+) concentration profiles in rabbit cells, in which the sarcoplasmic reticulum is disabled. Our model results agree with experimental measurements of global Ca(2+) transient in myocytes loaded with 50 μM Fluo-3. We found that local Ca(2+) concentrations within the cytosol and sub-sarcolemma, as well as the local trigger fluxes of Ca(2+) crossing the cell membrane, are sensitive to details of t-tubule micro-structure and membrane Ca(2+) flux distribution. The model additionally predicts that local Ca(2+) trigger fluxes are at least threefold to eightfold higher than the whole-cell Ca(2+) trigger flux. We found also that the activation of allosteric Ca(2+)-binding sites on the Na(+)/Ca(2+) exchanger could provide a mechanism for regulating global and local Ca(2+) trigger fluxes in vivo. Our studies indicate that improved structural and functional models could improve our understanding of the contributions of L-type and Na(+)/Ca(2+) exchanger fluxes to intracellular Ca(2+) dynamics

Warren McCulloch and the British cyberneticians

Warren McCulloch was a significant influence on a number of British cyberneticians, as some British pioneers in this area were on him. He interacted regularly with most of the main figures on the British cybernetics scene, forming close friendships and collaborations with several, as well as mentoring others. Many of these interactions stemmed from a 1949 visit to London during which he gave the opening talk at the inaugural meeting of the Ratio Club, a gathering of brilliant, mainly young, British scientists working in areas related to cybernetics. This paper traces some of these relationships and interaction

Molecular and Subcellular-Scale Modeling of Nucleotide Diffusion in the Cardiac Myofilament Lattice

AbstractContractile function of cardiac cells is driven by the sliding displacement of myofilaments powered by the cycling myosin crossbridges. Critical to this process is the availability of ATP, which myosin hydrolyzes during the cross-bridge cycle. The diffusion of adenine nucleotides through the myofilament lattice has been shown to be anisotropic, with slower radial diffusion perpendicular to the filament axis relative to parallel, and is attributed to the periodic hexagonal arrangement of the thin (actin) and thick (myosin) filaments. We investigated whether atomistic-resolution details of myofilament proteins can refine coarse-grain estimates of diffusional anisotropy for adenine nucleotides in the cardiac myofibril, using homogenization theory and atomistic thin filament models from the Protein Data Bank. Our results demonstrate considerable anisotropy in ATP and ADP diffusion constants that is consistent with experimental measurements and dependent on lattice spacing and myofilament overlap. A reaction-diffusion model of the half-sarcomere further suggests that diffusional anisotropy may lead to modest adenine nucleotide gradients in the myoplasm under physiological conditions